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Abstract: The measurements of the properties of the Higgs boson still leave room for 
a non minimal scalar sector. Extensions of the Standard Model typically involve multiple 
neutral Higgs helds which can interfere among themselves. We show that these interference 
effects can be substantial taking as example the one Higgs Singlet Model, the simplest 
renormalizable addition to the SM. 
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1 Introduction 

All mesurements performed in LHC Run I relative to the resonance discovered at about 
125 GeV [1, 2] are consistent with the hypothesis that the new particle is indeed the 
Standard Model Higgs boson. While the mass is already known to an astonishing three per 
mill from the latest published CMS measurement [3], the signal strengths /r* = cr*/cj^^, 
where i runs over the decay channnels, are known to about 10 to 20% [3-5]. This leaves 
room for modifications of the SM with a more complicated Higgs sector provided they are 
consistent within experimental errors with the data. A larger Higgs sector implies that 
additional scalar states are present in the spectrum. Direct searches have provided limits 
on the existence of new spin zero particles and on the strengths of their couplings [6]. 
The larger luminosity and energy in Run H will allow more precise measurements of the 
already discovered Higgs properties and extend the mass range in which other scalars can 
be searched for. 

The simplest renormalizable extension of the SM is the one Higgs Singlet Model 
(IHSM). It introduces one additional real scalar field which is a singlet under all SM 
gauge groups. The IHSM has been extensively investigated in the literature [7-28], how¬ 
ever, to our knowledge, no public MC implementation of the model is available. In this 
note we present such an implementation using FeynRules for the derivation of the vertices 
and Madgraph 5 for the generation of the amplitudes. We then discuss the simple case of 
Higgs production via gluon fusion at the LHC. 

2 The Singlet Extension of the Standard Model 

The singlet extension of the SM is defined by adding to the standard Lagrangian the 
following gauge invariant, renormalizable term: 

Jf, = di^Sdf,S - (2.1) 
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where S' is a real SU{2)l ® U{1)y singlet and <1>, is the SM Higgs weak isospin doublet. 
Here and in the following we adopt the notation of ref. [23]. A Z 2 symmetry which forbids 
additional terms in the potential is assumed. A detailed discussion of the IHSM without 
Z 2 symmetry can be found in ref. [26]. 

The neutral components of these fields can be expanded around their respective Vac¬ 
uum Expectation Values (VEVs) as follows: 


$ = 


Vd + f + iG^ 

V 71 


The minimum of the potential is achieved for 


S 


Vs + S° 

~vr 


( 2 . 2 ) 


A = Moi + ( 2 . 3 ) 

provided 

Ai, A2 > 0 ; 4A1A2 — A| > 0 . ( 2 . 4 ) 

The mass matrix in the gauge basis can be diagonalized into the (tree-level) mass basis 
introducing new fields: 

/i = cos a — sin CK and H = sina + cosa (2.5) 

with — I < a < |. 

The masses are 


m 


h,H 


= Ai + A 2 T |Ai — A 2 + tan^(2a), tan(2a) = 


hvdVs 


Aiu^ - A2 u2 ’ 


( 2 . 6 ) 


with the convention m'jj > m\. They correspond to a light [h] and a heavy [H] CP-even 
mass-eigenstate. 

The Higgs sector in this model is determined by five independent parameters, which 
can be chosen as 

nih, ruH, sin a, Vd, tan/3 = Vd/vs , (2.7) 

where the doublet VEV is fixed in terms of the Fermi constant through uj = Gj} 1^/2. 
Furthermore one of the Higgs masses is determined by the LHC measurement of 125.02 
GeV. Therefore, three parameters of the model are presently free. 

As only the doublet component, before mixing, can couple to the fermions (via ordinary 
Yukawa interactions) and the gauge bosons (via the gauge covariant derivative), all of the 
corresponding Higgs couplings are rescaled universally, yielding 


gxxs = gxfhi'^ + ^xs) with l-FAa;^ = 


cos a s = h 
sincK s = H 


( 2 . 8 ) 


- 2 - 









r(H^ hh) 


hh Width 


BR(H^ hh) 




Mh [GeV] Mh [GeV] 


Figure 1. On the left: the partial width for the H ^ hh process for Sq, = 0.3, 0.2, 0.1 and 
tan/3 = 1.0. On the right: the corresponding branching ratios. 


9xxs\S 2 — 9xxhhi^ “I" ^a;si)(l + ^XS2)- (2-9) 

where xx represents a pair of SM fermions or vectors. 

The couplings in the scalar sector involve tan /3 and are slightly more complicated. As 
an example we reproduce the triple scalar vertices in terms of the independent parameters 
mentioned above (sq, = sin a, Cq = cos a): 


^hhh 

3ieM^ o o 

(2.10) 

VhHH 

- tan/3}{Mf + 2Mjj} 

(2.11) 

yHhh 

- +«« tanmmi + 

(2.12) 

Vhhh 
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(2.13) 
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The tree level partial width for the decay of the heavy scalar into two light ones reads: 


r{H hh) 


e^Mfj f 4M2\ 

1287rV J 
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SaCo, 


(Ca + Sa tan /3) 


2 


(2.14) 


In fig. 1 we show the partial width and the corresponding branching ratios as a function 
of the heavy Higgs mass for Sa = 0.3, 0.2, 0.1 and tan (3 = 1.0. The BR is computed as the 
ratio of the lowest order width in eq.(2.14) to the total width given in ref. [29] multiplied by 
Sq. The BR rises sharply above the kinematical threshold and, for the parameter range we 
have considered, remains larger than 25% up to Mjj = 1 TeV. This raises the tantalizing 
prospect of a relatively abundant production of heavy Higgses followed by their decay into 
two light ones if Mjj > 2Mh- See ref. [26] for a detailed discussion in the framework of the 
ISHM. 


3 Limits on the parameters 

The strongest limits on the parameters of the IHSM come from measurements of the 
coupling strengths of the light Higgs [3-6], which dominate for small masses of the heavy 
Higgs, and from the contribution of higher order corrections to precision measurements, in 
particular to the mass of the W boson [24] , which provides the tightest constraint for large 
Mh- The most precise result for the overall coupling strength of the Higgs boson from 
CMS [3] reads 


fi = a/asM = 1-00 ±0.13. (3.1) 

Therefore the absolute value of sin a cannot be larger than about 0.4. This is in agreement 
with the limits obtained in ref. [24, 27] which conclude that the largest possible value for 
the absolute value of sin a is 0.46 for Mh between 160 and 180 GeV. This limit becomes 
slowly more stringent for increasing heavy Higgs masses reaching about 0.2 at Mh = 700 
GeV. 


4 Interference effect and simulation tools 


The focal point of this note is the interference between the two Higgs fields. In general, 
any amplitude involving a single Higgs exchange can be written as 


4 = 4' 


+ 


- m 2 + iThMh q^-M%± iThMh 


+ 4o = 4i + 4o 


(4.1) 


where Aq does not involve the scalar fields. The real parts of the two propagators in- 
tefere destructively for M| < < M‘^ and constructively for < M| and M^ < q'^. 
This phenomenon has already been noticed in ref. [25] where however it was dismissed as 
numerically irrelevant. 

As will be argued in sect. 5, the interference effects can be substantial, their relevance 
increasing with Mh- The shape of the the heavier Higgs peak is also strongly affected and 
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this will need to be taken into account in any search for additional scalars and eventually 
in the measurement of their properties. 

These features are neglected by any prediction based on the narrow width approxima¬ 
tion, or equivalently on a production times decay approach. 

The relevance of the interference term does depend on the relative size of 
in eq.(4.1) and is expected to be significant for processes in which the resonant part of the 
amplitude is large as in gluon fusion for large values of the heavy Higgs mass. 

Clearly, this interference between different Higgs fields is not a peculiarity of the Singlet 
Model. It will indeed occur in any theory with multiple scalars which couple to the same 
set of elementary particles, albeit possibly with different strengths. 

In order to allow for the Monte Carlo simulation of the IHSM we have used Feyn- 
Rules [30, 31] to prepare a UFO file [32] for the model, which can be imported, as we 
did, in MadGraph 5 [33] and many other general purpose MC tools. It enables the sim¬ 
ulation at tree level of any process in the IHSM. The UFO file can be downloaded from 
http://personalpages.to.infn.it/~maina/Singlet. 

The gluon fusion channel in MadGraph requires particular care. MadGraph treats the 
gluon-Higgs effective vertex in the narrow width approximation, through an expansion of 
the top loop amplitude in powers of M|/M|,p, which is unsuitable in the present context 
and which fails for Higgs masses above the t-tbar threshold. The appropriate effective 
vertex Vggh must be introduced by hand in the matrix element: 


with 



as{V2GF)"^Ml^ 

TTS 


(1 




C{t) 


— 2 arcsin(l/-y/r)^ t>1 

5 (log(l^) -ivr) r<l 


(4.2) 


(4.3) 


where r = AM^^/s, (3 = y/l — r and s is the square of the sum of the momenta carried by 
the gluons which in general is not equal to the Higgs mass squared. 

The tensor structure of the vertex is already taken care of by MadGraph. 


5 Results 

As an example we have studied Higgs production through gluon fusion in the four electron 
channel, gg ^ h, H ^ ZZ —?• 4e, at the LHC with a center of mass energy of 13 TeV. 

We neglect the non resonant contribution given by gg —>■ ZZ and hy gg ^ ZZ through 
a quark box amplitude. This continuum is known to be large and in particular there is 
a non negligible interference between the box contribution and Higgs production through 
the heavy quark three point loop. These terms are well known [34-47] and are essential 
for accurate phenomenological predictions. The region of large invariant masses of the 
four final state leptons in gg —)• ZZ, WW —)• Al has been studied in detail in ref. [47] 
in the SM with a Higgs of 125 GeV. Above the light Higgs peak, the differential cross 
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section is dominated by the continuum for M/^i < 2Mtop- The box contribution however 
drops more rapidly with increasing M 41 than the Higgs mediated one. Unitarity requires 
the interference between these two components to be negative and while it is essentially 
negligible below the top threshold it becomes more relevant for larger masses and exceeds 
50% of the incoherent sum of the two contributions in the one TeV range. 


41 invariant mass 



Figure 2. gg ^ h,H ^ ZZ — >■ 4e at the LHC with a center of mass energy of 13 TeV. Upper 
part:four lepton invariant mass distribution for = 125 GeV, Mh = 400 GeV and Sa = 0.2. 
The red line is the full result (u). The violet histograms {(Jh) shows the SM prediction with Higgs 
couplings scaled by cos a. The blue line {(Jh) gives the result when the light Higgs diagrams are 
neglected while the H ^ hh contribution to Th is retained. The black histogram (ctat/) refers to 
the sum of the violet and blue lines and corresponds to neglecting the interference effects. Lower 
part: the ratio . 
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Figure 3. Four lepton invariant mass distribution for Mh = 125 GeV, Mh = 400 GeV, Sa = 0.3. 
The meaning of the various histograms is as in fig. 2. 

Our main results are shown in figs. 2-5 and tab. 1. No cut has been applied to the final 
state. We have used CTEQ6L1 parton distribution functions [48]. The ratio of vacuum 
expectation values tan j3 has been taken equal to one. 

The width of the light and heavy scalar are 

T;, = Th = V^^{Mh)sI + V{H ^ hh). (5.1) 

This corresponds, using eq.(2.14) for the H ^ hh width and ref. [29] for the SM Higgs 
width, to Til/ = 1.77(4.08) GeV for Mh = 400 GeV, Sq, = 0.2(0.3); Th = 15.80 GeV for 
Mh = 600 GeV and Sq = 0.3; Th = 16.69 GeV for Mh = 800 GeV and Sq, = 0.2. We have 
assumed Mh = 125 GeV which corresponds to r'^^(125) = 4.03 MeV. 
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Figure 4. Four lepton invariant mass distribution for Mh = 125 GeV, Mh = 600 GeV, Sa = 0.3. 
The meaning of the various histograms is as in fig. 2. 



200 GeV < Mil < 1 TeV 

M}{—2h GeV < Mu < Mh+25 GeV 


a 

<7/i + CTh 

cth 

a 

O'h + cth 

CTH 

Mh= 400 GeV, Sa=0.2 

72.95 

70.96 

26.00 

32.09 

32.13 

25.55 

Mh= 400 GeV, Sa=0.3 

101.48 

96.51 

55.99 

60.36 

59.70 

53.74 

Mh=600 GeV, Sa=0.3 

48.44 

52.52 

11.99 

11.51 

11.96 

9.97 

Mh=800 GeV, Sa=0.2 

43.96 

45.96 

1.00 

1.57 

1.46 

0.65 


Table 1. Cross sections in ab for gg ^ h, H ^ ZZ —>■ 4e at the LHC with a center of mass energy 
of 13 TeV. 
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Figure 5. Four lepton invariant mass distribution for Mh = 125 GeV, Mh = 800 GeV, Sa = 0.2. 
The meaning of the various histograms is as in fig. 2. 

We show the region around the heavy Higgs peak where the interference affects are 
most prominent. The invariant mass distribution in the neighborhood of the light Higgs 
resonance is unaffected within the accuracy of the present simulation. In all figures the 
violet histogram is the result obtained including only the light Higgs with SM couplings 
scaled by which we denote as ah- Since T/j is small ah ~ crsM{Mh)c^ for large M^i. 
The blue line shows the cross section when only the heavy Higgs is present, which will be 
referred to as an- The H ^ hh contribution to T/f is retained and therefore an is not 
equal to asM{^H)Sa- The full result, a, is shown in red in fig. 2 for Sa = 0.2, Mh = 400 
GeV, in figs. 3, 4 for Sa = 0.3 with Mh = 400 GeV and Mh = 600 GeV respectively and 
finally in fig. 5 for Sa = 0.2, Mh = 800 GeV. The black histogram displays the sum of 
the blue and violet lines, a^i = (^h + and corresponds to neglecting the interference 
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between the two scalars. The fractional size of the correction to ajsji is displayed in the 
bottom part of the figures where the ratio of the interference term and is 

shown. 

Contrary to naive expectations the light Higgs contribution is non negligible outside 
the peak region [47] and the interference effect is substantial. While the details depend 
obviously on the mass of the heavy Higgs and on the mixing angle, we find a decrease of 
10 to 20% of the differential cross section at invariant masses of the four leptons of 300 
GeV for Mh = 400 and of about 20% for Mu = 350, Mh = 600 and for M^i = 500, 
Mh = 800. This depletion becomes more pronounced as increases and reaches a dip 
which is almost two orders of magnitude smaller than the predictions which neglect the 
interplay of the two Higgs fields. The interference stays negative for four lepton masses 
below the heavy Higgs mass and then turns positive. It attains a maximun in which the 
true value is about a factor of two larger than aMi and then slowly decreases. The position 
of the peak is shifted to slightly larger masses. At four lepton masses about 200 GeV larger 
than Mh the interference still amounts to about 20 to 40% of a mi- Because of unitarity, 
for very large M^i the full cross section a must approach the SM result (Tsm- In this region, 
where the width of the two Higgses can be neglected, asM = Ic% = ^h! 

tab. 1 shows the cross section in ab for two mass intervals: 200 GeV < M 41 < 1 TeV, 
which roughly coincides with the range employed so far by the experimental collaborations 
to set limits on the presence and couplings of additional scalars, and M//—25 GeV < M^i < 
Mh+‘^5 GeV, as an indication of the possible effects on an analysis in smaller mass bins 
which requires high luminosity. In the hrst case, the contribution of the heavy Higgs is a 
relatively small fraction of the Higgs production cross section in gluon fusion. Furthermore, 
the interference effects in the 200 GeV < M 41 < 1 TeV depend crucially on the heavy Higgs 
mass. For Mh = 400 the exact result is larger than the incoherent sum + an- The 
long tail for M/^i > Mh gives a larger contribution than the intermediate region 200 GeV 
< M 4 / < Mh- On the contrary, for larger heavy Higgs masses we have a < ah + an- 
The negative interference in the intermediate region outweights the positive contribution 
at larger masses. In the smaller range, Mh—25 GeV < M^i < Mh+25 GeV, the positive 
and negative contributions are very close for Mh = 400 GeV and the full result is in 
rough agreement with ah + an- For this value of Mh and mass interval, the heavy Higgs 
contribution is significantly larger than the light Higgs one. For Sa = 0.3, Mh = 600 GeV 
the exact cross section is about 4% smaller than the incoherent sum while for Sq = 0.2, 
Mh = 800 GeV it is approximately 8% larger. We see that for large values of the heavy 
Higgs mass the interference effects are non negligible even on a restricted mass range. 

Our results have no pretense to be a complete prediction. They need to be validated 
with the inclusion of the continuum contribution and of higher order corrections. It should 
be noticed that, for Higgs decay to color neutral hnal states, all relevant amplitudes in QGD 
will have the structure of eq.(4.1) and therefore interference effects between the scalar fields 
will not be spoiled by QGD corrections. 

The interference with the quark box amplitude deserves more care. The gluon-gluon 
continuum term does not involve scalar exchanges and therefore it cannot be cast in the 
form of eq.(4.1) and could in principle dilute the effect. However, the interference between 
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the box diagrams and the light Higgs mediated ones is always negative, while the ampli¬ 
tude with a heavy Higgs exchange changes sign at the resonance. As a consequence the 
continuum and the heavy Higgs term are expected to be in phase for M/j < M 4 ; < Mh 
and out of phase for > Mh- A detailed study of this topic is in preparation. 

6 Conclusions 

In any theory with multiple neutral Higgs which couple to the same set of elementary 
particles the scalars are expected to interfere. We have shown in the case of Higgs produc¬ 
tion through gluon fusion in the ISHM that the interference effects can be significant and 
cannot be neglected when aiming for high accuracy predictions. 
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